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The present invention relates to a method of velocity filtering seismic signals, and to an installation for imple- 

The invention relates to studying underground formations, in particular by the so-called "Vertical Seismic 
Profile" technique in which a seismic wave detector is placed in a borehole successively at different depths, 
with seismic waves being emitted from a source on the surface and with the signals produced by the detector 
being recorded. 

The essential purpose of these measurements is to determine the reflective horizons or reflectors which 
are situated deeper than the bottom of the borehole by analyzing the waves reflected on these reflectors and 
rising towards the detector, which waves are called "upgoing" waves. 

However, the detected waves comprise not only the upgoing waves, but also waves known as "downgoing" 
waves which are waves that have propagated directly from the source to the detector, waves which have been 
subjected to multiple reflections, and interfering waves of various kinds. 

By putting together all of the recorded signals in a single document, it is possible to detect coherencies 
between the various traces, but given the multiplicity of superposed components in each signal, such a docu- 




i a set of signals r 

in 1 and n) is described in patent document FR-A-2 494 450. This filtering cc 
the m (m < n) first signals g,(t) recorded at levels z, by a time t, where 2 < i< m relative I 
level Z! in order to align the downgoing wave components and produce a first signal by adding together the 
signals shifted in this way. Similarly, a second signal is produced by adding together recorded signals g,(t) for 
1< i < m after delaying the signals g,(t) for 2 < i < m by a time t, relative to the signal g,(t). The first and second 
signals are then combined in order to produce a signal u,* representative of an optimum estimation of the upgo- 
ing wave component. Signals u 2 * u k * are produced in the same mannerfrom setof m recorded signals {g 2 (t), 

...,9n. + 1«t)} {SkTO g m *i,-i<t)>..... 

This method suffers from the drawback of reinforcing, in practice, only those waves which are aligned in 
the recorded signals after they have been shifted, whereas waves having different velocities simply do not show 
up. 

There also exists the so-called F-K velocity filter method which serves to separate upgoing waves from 
downgoing waves. Reference can be made to "Vertical Seismic Profiling" by B.A. Hardage, Geophysical Press, 
1983, pp. 175-179 for a detailed description of this method. 

In this method, the set of signals g,(t) recorded by the detector at depths z, for 1 < i < n Is taken to be a two 
dimensional signal g(z, t) and the following successive operations are performed : 

a) The two-dimensional Fourier transform G(k, w) of the signal g(z, t) is calculated where k and w are the 
dual variables of z and t respectively in the Fourier transform ; the signal G(k, w) has the advantage that 
the upgoing waves are to be found in the quadrant k > 0, w > 0 and the downgoing waves are in the quad- 
rant k < 0, w > 0 ; 

b) The coefficients contained in the quadrant k < 0, w > 0 are set to zero or are multiplied by a small factor, 
e.g. 10-3; and 

c) The inverse two-dimensional Fourier transform is calculated and in the resulting signal g u (z, t) the 
downgoing waves are highly attenuated. 

This method suffers from problems due to side-effects because of the sharp cut-off between those coeffi- 
cients which are modified in operation b) and the other coefficients which are not modified. In particular, this 
gives rise to artifacts and to spectrum folding during operation c). 

Further, the Fourier transform using each of the variables z and t requires a certain number of samples to 
be available. This is no problem for the variable t since there are generally several thousand points available 
along the t-axis, given that each signal is recorded'over a period of several seconds and that the sampling period 
is about 1 ms. In contrast, the number of points available along the z-axis depends on the number of levels at 
which signals are recorded, and this is much less. 

in orderto perform the Fouriertransform on the variablez, it is necessary to have at least 64 different record- 
ing levels, with the gap between consecutive levels being no greater than a value related to the highest fre- 
quency in the received signal and to the lowest propagation velocity, with a typical value being about 10 meters 
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measurement costs. However, this method suffers from the same defects of artifacts and spectrum folding as 
the F-K method. 

id of spectrum folding which appear, in 
to make it possible to perform 
velocity filtering on a small number of recorded signals in order to reduce the cost of on-site measuring and to 



More precisely, the present invention provides a method of seismic exploration by filtering a two-dimen- 
sional signal g(z, t) built up from a set of signals g,(t) for 1 < i< n as produced by at least one seismic wave 
detector placed at different depths z, z z„ in a borehole, each of said signals being produced in response 



A = 1/2 fjd + e.B.H,.DJ 



B is a normalization factor whose value depends on the signal to which the operator A is applied ; 
H t is the one-dimensional Hilbert operator with respect to the variable t ; 
D is a differentiation operator with respect to the variable z ; and 
e is an integer equal to + 1 or - 1 ; 

al signal g(z, t), with the resulting signal being 
sd by upgoing waves when e = + 1 and by downgoing waves when e = - 1 . 
The normalization factor B may be defined, for a function g, as follows : 



J IIVDf 



where || || indicates the L2 norm of the function. 

It will thus be understood that when the operator A Is applied recursively to a signal, the normalization factor 
B must be recalculated for each iteration as a function of the signal to which the operator A is applied. 

The Hilbert operator H, may be defined as follows : 

H l (g(t» = FT-i0.sgn( 0) ).FT(g(t))] 



g(t) is a function of the variable t ; 
FT-' is the inverse Fourier transform ; 
I is the complex operator such that j* = - 1 
sgn is the sign function which is equal to + 




the number of samples available is commonly several th 
to the variable z does not require a large number of recordings z, to be av 
or D z applied to a signal g(z) is defined by : 



D z [g(z)] -]jH 



g(z+Zm-l) - g(z-2m+l) 



waves having velocities lying between two predetermined velocities Vc, and Vc* This is done as follows : the 
recorded signals a(t), 2 < i < n constituting the two-dimensional signal g(z, t) are shifted in time relative to one 
another so that the components of waves of velocity Vc, appear at the same instant in each signal ; the shifted 
two-dimensional signal is then filtered in accordance with the invention and the resulting signal is time shifted 
to cancel the initial shift. The same operation is performed for the two-dimensional signal g(z, t) for the velocity 
VC2 in order to produce a second two-dimensional signal. The difference between the two two-dimensional sig- 

is whose velocities 




I a vertical seismic pro- 
file ; 

Figure 2 shows in highly simplified manner the paths of upgoing and down seismic waves relating to two 



Figure 3 shows the positions of the upgoing waves U and the downgoing waves D in the o-k plot of the 
F-K transform domain ; 

Figure 4 is a diagram of the transfer function of the operator R in the oo-k plot of the F-K transform domain; 
Figure 5 is a diagram of the transfer function of the operator A in the a-k plot of the F-K transform domain; 
Figures 6 to 9 are respectively a signal recorded by means of an installation in accordance with Figure 1 ; 
and the upgoing wave signal after one (Figure 7), two (Figure 8) and three (Figure 9) applications of the 
operatorAwithl = 1 ; 

Figures 1 0 to 1 2 respectively show the signal recorded by i 
ure 1 , and an upgoing wave signal and a downgoing i 
Figure 13 is a flow chart of processing means in an i 

Figure 14 is a plot of the co-k transform showing the zone cc 
between the velocities Vc, and Vcj ; and 
Figures 15a to 15e are plots showing the filtering of a detected signal g(z, t) for eliminating waves having 
a velocity less than Vc, from this signal. 

Figure 1 shows an installation for seismic exploration comprising a drilling derrick 2 disposed over a 

A downhole apparatus 8 including a seismic wave sensor such as a geophone is suspended down the 
borehole 4 by means of an electric cable 10 which passes over sheaves mounted on the drilling derrick. The 
cable 10 is used for displacing the apparatus 8 within the borehole and simultaneously serves to transmit the 
detection signals produced by the sensor towards surface equipment 12. In conventional manner, the surface 
equipment 12 comprises a winch for winding the cable and means for determining the depth at which the dow- 
nhole apparatus 8 is located, together with means for processing and recording the detection signals transmit- 
ted via the cable 10. 

The downhole apparatus 8 may be conventionally constituted by a moving anchor member 14 capable of 
anchoring in the wall in order to ensure appropriate contact between the sensor and the wall of the borehole. 
This takes place only when the apparatus 8 has arrived at a depth at which a measurement is to be performed, 
with the surface equipment 12 then transmitting a message in order to cause the element 14 to move into the 
anchoring position. 

A source of seismic waves 16 is located on the surface at a given distance from the head of the borehole. 
This source may be constituted by any suitable device such as an air gun. The source 16 is controlled from 
the surface equipment 12 in order to fire a shot when the downhole apparatus 14 is anchored at one of the 



Several shots are generally fired for each depth, and all of the signals produced in this way are stored and 
added together in the surface equipment 12 so as to obtain a signal for each level in which the significant com- 
ponents are reinforced to the detriment of random noise components. This signal is referred to below as the 
detection signal. 

The set of detection signals g,(t) obtaii 
recordings (known as seismic traces) whic 

Figure 2 shows highly simplified versions of the paths of thi 
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different levels. Reference 1.8 designates a reflecting horizon or -reflector" which is deeper than the bottom 20 
of the borehole. The reflector 18 is constituted by the interface between two layers of significantly different 
acoustic impedances. , „ 

The sensor placed at A, at level z, receives a sound wave coming directly from the source and called the 
downgoing wave, and also receives a wave reflected by the reflector 18. and called the upgoing wave The 
upgoing wave needs to be shown up in order to be able to determine the depth of the reflector and also its 
reflection coefficient. 

The detection signal g,(t) obtained at level z, is thus the sum of a downgoing wave component and an upgo- 
ing wave component, and in addition it includes noise components, some of which come from multiple reflec- 
tions on reflectors situated between the surface and the bottom of the borehole. 

The sensor placed at A k at level z k similarly receives a downgoing wave and an upgoing wave, however 



. .ie method of the invention consists in processing the two-dimensional signal g(z, t) in order to isolate 
eitherthe upgoing waves, or else the downgoing waves. „„„«,„,,„ t „ , ha 

The Invention lies firstly in the definition of two functions g u and g d which correspond respectively to the 
upgoing waves and to the downgoing waves, and secondly in the use of approximate expressions for these 
functions in order to provide a velocity filter. 
We begin by showing the functions : 

g„ = «2[g + H2(g)]and 
g d = 1/2fe-H2(g)] 

where H2 is the two-dimensional Hilbert transfo.ro respectively containing the upgoing waves and the downgo- 
ing waves of the recorded signal g(z, t). 

This is achieved by using the known properties of the Fourier transform. It is known that in the F- 
domain, the upgoing waves lie in the quadrants for which the variables m and k have opposite signs, wl 



In ^Figure 3 'plot the upgoing waves U lie in quadrants II and IV, and the downgoing waves D lie in quad- 
;s I and III. 
Using the d< 

FTTH(f)]M = j.sgn( M ).FT[f](o>) (1) 



rants I and lit rf ^ Fourferand ofthe H i|bert transforms, the following equation can be deduced 



FT and H are the one-dimensional Fourier and Hilbert transforms, respectively : 
f is a function; 

j is the complex operator such that f = - 1 ; and 
sgn is the sign function which is defined by 

sgn( ffl )-+1 ifto>0 
sgn(<o) = -1 ifo><0 
sgn(d>) = 0if(» = 0. 

The next relationship is easily derived from equation (1) : 

FT2[H2(g)K<D. k) = - sgn(co, k).FT2[g](<i>, k) 

where FT2 and H2 are the two-dimensional Fourier and Hilbert transforms respective 
FT2[gJ(o>, k) can then be expressed in the following form : 



FT2[g u ](<D, k) = 1/2 (FT2[g](<», k) + FT2[H2(g)](<o, k)) 
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The coefficient (1 - sgn(o, k))/2 is equal to 1 in quadrants II and IV and is equal to 0 in quadrants I and III 
of the Figure 3 plot, thus showing that the function g u contains upgoing waves only. 
It can similarly be shown that 



which means that the function g d is equal to the signal g in quadrants I and III, and is equal to 0 in quadrants 
II and IV, i.e. that the function g t contains downgoing waves only. 

Velocity filtering in accordance with the invention thus consists in applying the operator R equal to 1/2Jld 
+ sH2], to the recorded signal g(z, t), with Id being the identity operator and s = - 1 or + 1 depending on whether 
20 the upgoing waves or the downgoing waves are to be obtained. 

The two-dimensional Hilbert transform H2 comprises a transform H t with respect to the variable t and a 
transform H z with respect to the variable z. In order to apply these transforms directly, it is necessary for the 
number of samples in the signal along the variable in question to be large, for example not less than 100. 

There is no problem for the H, transform since several thousand points are generally available along the 

In contrast, the number of points along the z-axis, i.e. the number of recorded signals gj(t) with 1 < i < n is 
limited for reasons of cost. The operator H z cannot therefore be applied directly. 

In accordance with the invention, the operator H z is replaced by a differentiation operator D z , i.e. an operator 
operating on a limited number of points along the z-axis and which constitutes an approximation to the operator 
30 H z , and the operator R modified in this way is applied recursively. 

More precisely, the operator R equal to 

1/2[ld + eH,HJ 

35 is replaced by the operator A which is equal to 

1/2[ld+eB.H,.D z J 

where B is a normalization coefficient whose value depends on the signal g to which the operator A is applied 
40 in such a manner that the norm of the signal A[g] is equal to the norm of the signal g. 
Particular implementations of the operators H, and D z are now described. 

The Hilbert transform H, may be implemented by any of the methods known to the person skilled in the art. 
In particular, the conventional formulation of the Hilbert transform using a convolution may be used : 

45 H,(g) = g.CONV.(1/t) 

However, it is preferable to use the following equation deduced from equation (1) : 

H l (g) = IFT,D.sign( 0 >).FT,(g)] 

where IFT, is the inverse Fourier transform. 

In this equation, the Hilbert transform of the function £ is obtained by calculating the Fourier transform of 
the function g_, multiplying the result by jisgnfro), and calculating the inverse transform. This method has the 
advantages of being very fast and of providing a very accurate result. As for the operator D„ it is advantageous 
55 to use an expression of the following form : 



6 



g(z+2m-l) - g(z-2m + l) 



where I is an integer whose value is not more than a few units. Fori = 1, D cc 
for I = 2 it constitutes a seven-level filter. 
10 It has been observed experimentally that filtering is already very effective for I = 1 . 

In this particular case, the filtering process of the invention may be explained, from a theoretical point of 

From equation (1) it can be deduced that the Hilbert transform H, has j.sgn(m) as its transfer function. Simi- 
larly, H 2 has j.sgn(k) as its transfer function. From this it can be deduced that the operator R has a transfer 
is function T(R) equal to : 



30 Figures 4 and 5 are diagrams of the transfer functions T(R) and T(A) in the co-k plot of the F-K transform 

domain. The coefficients shown in the quadrants are the multiplicative coefficients applied to each signal 
sample In the quadrant 

In Figure 4, the value of the multiplicative coefficient switches abruptly from one value to another when 
moving from one of the quadrants into an adjacent quadrant. As mentioned above, this gives rise to defects in 
35 the filtered signal, and in particular it gives rise to spectrum folding. 

In contrast, it can be seen that the value of the multiplicative coefficients shown in Figure 5 changes con- 
tinuously on moving from one quadrant to another. This explains why a signal filtered in accordance with the 
method of the present invention does not suffer from spectrum folding. 

It is also recalled that in accordance with the invention the operator A is applied recursively to the recorded 



40 signal g(z,t). Thus, in fact, it is the operator Ap for integer pth 

are the coefficients shown in Figure 5 raised to the power p. It is clear that A? gets closer to R for increasing 
E- 

It has been observed experimentally that, for upgoing waves, the filtered signal is already usable with p = 
1, and becomes very good with p = 2. As for downgoing waves, the filtered signal is already good with p = 1 
45 since the signal/noise ratio for these waves is h igher than for upgoing waves. 

When filtering to show up the upgoing waves, it is possible to apply bandpass filtering over the band [2 Hz, 
50 Hz] after each application of the operator A In order to remove residual high frequency downgoing waves. 

Figure 6 shows a two-dimensional signal g(z, t) recorded by means of an installation such as that shown 
in Figure 1. The upgoing waves are not apparent in this raw signal. 
so Figure 7 shows a signal obtained after applying three-level filtering to the Figure 6 signal with an operator 
A for which e = + 1 and I = 1. The upgoing waves show up much more clearly, while the downgoing waves are 
attenuated. 

The signal shown in Figure 8 is the result of reapplying the operator A on the Figure 7 signal. The upgoing 
waves are nearly perfect and the residual downgoing waves can hardly be seen. 
55 Figure g shows the result of applying the operator A to the Figure 8 signal. Figure 9 is pratically identical 
to Figure 8 which means that filtering convergence has already been obtained in Figure 8. 

Figure 10 shows another signal recorded by means of an installation such as that shown in Figure 1 . This 
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with an operator A in which e =- 1 and I = 1. The corresponding upgoing and downgoing signals are shown 
respectively in Figures 11 and 12. 

Figure 13 is a flow chart applicable to an installation in accordance with the invention for filtering a two- 
dimensional signal g(z, t). 

This installation comprises a memory 22 for receiving the detection signals g,(t), for 1< i < n constituting 
the two-dimensional signal g(z, t). 

This memory is controlled to provide the necessary signals g,(t) to calculating means 24 in orderto produce 
the result of applying the operator D z to the signal g(z, t). E.g., when the operator D z is defined by DJIifz)] - 
h(z - 1) - h(z + 1) for some function h, the memory is controlled to provide the calculating means 24 with pairs 
of values (g,(t), g l+2 (t)) for 1< i < n - 2. 

The signal DJg(z, t)] delivered by the calculating means 24 is then received by calculating means 26 wnicn 
produces the Hilbert transform relative to the variable t This Hilbert transform may be performed in any known 
manner, and in particular, as mentioned above, it may be performed using the Fourier transform in the following 

H,(g) = IFT 1 D.sgn(< i >).FT,(g)] 

The signal H,.D 2 .[g(z, t)] delivered by the calculating means 26 Is applied to one of the inputs of calculating 
means 28 ; the signal contained in the memory 22 is applied to the other input of the calculating means.28. 
This calculating means 28 calculates the value of the coefficient B which appears In the operator A of methods 
in accordance with the invention to make the operator A conserve the amplitude of the signal applied thereto. 

The coefficient B delivered by the calculating means 28 is then multiplied by the signal delivered by the 
calculating means 26 in a multiplier 30. The multiplier 30 also receives the coefficient e whose value is fixed 
by the user : e = + 1 if the upgoing waves in the signal g(z, t) are to be retained, or else c = - 1 if the downgoing 
waves are to be retained. 

The signal e.B.H t .D 2 .[g(z, t)] delivered by the multiplier 30 is received on one of the inputs to an adder 32 
whose other input receives the signal contained in the memory 22. The resulting signal is then applied to one 
of the inputs of a divide-by-two circuit 34. 

The means 24 to 34 constitute a particular embodiment for implementing the operator A In accordance 
with the invention, this operator is applied recursively. The number of times I that this operator is applied to the 
two-dimensional signal may be loaded, for example, into a counter 36 which controls a multiplexer 38. 

The input of the multiplexer is connected to the output from the divider 34 ; one of its outputs is connected 
to the memory 40 and its other output is connected to the memory 22. The counter is initially loaded with the 
value I and is decremented each time the operator A is applied. So long as the contents of the counter is not 
zero, the divider 34 is connected to the memory 22. When the counter reaches the value 0, the divider 34 is 
connected to the memory 40. 

A recorder 42 is connected to the memory 40 in order to produce a graphical recording of the processed 
two-dimensional signal as shown in Figures 7, 8, 9, 1 1, and 12. 

Figures 7, 8, 9, and 11 show signals which contain only upgoing wave components. In these signals, all 
of the upgoing waves are present regardless of their velocities. Similarly, the signal shown in Figure 12 includes 
all of the downgoing wave components, regardless of their velocities. 

It is possible to perform filtering over a fan so as to retain only those waves whose velocities lie in a pre- 
determined velocity band. Such filtering serves to eliminate noise components and components which are not 
applicable to the intended analysis to be performed on the detected signal g(z, t). 

For example, instead of conserving all of the upgoing w ■ - """•"-« 

plot, it may be desirable to conserve only those waves whc 
unshaded portion of the plot shown in Figure 14. 

Such fan filtering may be obtained in three stages by using the method of separating upgoing and downgo- 
ing components in accordance with the invention. The first stage consists in processing the signal g(z, t) to 
i produce a first signal g,(z, t) from which all waves having a velocity of less than Vc, have been eliminated. 

Similarly, the second stage consists in processing the signal g(z, t) in order to produce a second signal 
g 2 (z t) from which all waves having a velocity less than Vc z have been eliminated. 

Finally, the third stage consists in subtracting the signal g 2 (z, t) from the signal g,(z, t) in orderto produce 
a signal g 12 (z, t) which contains only those waves whose velocities lie between Vc, and Vcj. 

The first stage of such range filtering is described with reference to Figures 1 5a to 1 5e. In Figure 1 5a, three 
upgoing waves having different velocities Vc,, Vc,, and Vcj such that Vc > Vc, > Vcj are shown diagrammati- 

To begin with, the signals g,(t) for 2 < i < n constituting the two-dimensional signal g(z, t) are time shifted 



such that in each signal g,(t) for 1 < i < n the components of waves at velocity Vc, are aligned. As a result the 
upgoing wave of velocity Vc,, becomes a downgoing wave (Figure 15b). Such shifting is described in above- 
mentioned patent document FR-A-2 494 450, but for use in a different application. 

The velocity filtering method in accordance with the invention is then applied in ordertoshow up the upgoing 
wave components in the shifted signal. In the transform domain plot shown in Figure 15c, this gives rise to a 
high degree of attenuation for the waves contained in quadrants I and III. In the time domain plot of Figure 15d, 
waves of velocities less than Vc, are emphasized whereas waves of velocities greater than Vc, are highly 



The first stage terminates by shifting the signals g,(t) for 2 < i < n back again in order to cancel the time 
10 shift performed initially. This gives rise to the signal g,(z, t) which is shown in the Figure 1 5e plot 

The second stage is identical to the first stage except that the signals g,(t) for 2 < i < n are time shifted so 
that wave components of velocity Vcj are aligned in the signals g,(t) for 1< i < n. 

The time shifts performed at the beginning of the first and second stages are delays when the velocities 
Vc, and Vc 2 are positive (upgoing waves) ; the shifts would be advances If the velocities were negative (downgo- 



A signal g 12 (z, t) which predominantely contains waves whose velocities lie in the band ryc„ VcJ is obtained 
■ ■ the difference between the signal g,(z, t) produced in the first stage and the signal gj(z, t) 



1 . A method of seismic exploration by filtering a two-dimensional signal g(z, t) built up from a set of signals 
g,(t) for 1 < I < n as produced by at least one seismic wave detector placed at different depths z, Z| z„ 

as In a borehole, each of said signals being produced in response to seismicwaves being emitted, and said signals 
including upgoing wave components and downgoing wave components, said filtering consisting in separating 
said components and the method being characterized in that the operator A is applied to the two-dimensional 
signal g(z, t), where : 

30 A=1/2pd-fs.B.H t .D 2 ] 

in which : 

Id is the identity operator; 

B is a normalization factor whose value depends on the signal to which the operator A is applied ; 
35 H, is the one-dimensional Hilbert operator with respect to the variable t ; 
D is a differentiation operator with respect to the variable z ; and 
s is an integer equal to + 1 or - 1 ; 

said operatorbeing applied recursively to said two-dimensional signal g(z, t), with the resulting signal being 
predominantly constituted by upgoing waves when s = + 1 and by downgoing waves when s = - 1 . 
40 2. A method according to claim 1 , characterized in that the operator D z is defined by : 



- ■ L 



g(z + 2m-l) - g(z-2m+l) 



where g is some function and I is an integer. 

3. A method according to claim 2, characterized In that I is equal to 1. 

4. A method according to claim 2, characterized in that I is equal to 2. 

5. A method according to any one of claims 1 to 4, characterized in that the operator H, is defined by : 

H 1 fe(t)] = IFTrj.sgn(<D).FT(g(t)] 



g is soma function; 

FT represents the Fourier transform ; 
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sgn(a>) designates the sign of the variable to, which is the dual of t ; and 
] is the complex operator such that p = - 1. 

6. A method according to any one of claims 1 to 5, characterized in that the operator A is applied once only 
to the signal g(z,t). 

7. A method of velocity filtering a two-dimensional signal g(z, t) built up from a set of signals g,(t) for 1< i 

< n, said signals being produced by sound wave detectors placed at different depths z, z z„ down a 

borehole, each signal being produced in response to sound waves being emitted from a source on the surface, 
the velocity filtering serving to reinforce waves in the signal g(z, t) having velocities lying between predetermined 
velocities Vc, and Vc* said method being characterized in that it comprises the following operations : 

A) a) producing a first shifted two-dimensional signal by applying a time shift to the signal g,(t) for 2 < i < n 
in order to time align the wave components of velocity Vc, in the signals g,(t) for 1< i< n ; 

b) filtering said shifted signal using the method according to any one of claims 1 to 6 ; 

c) producing afirst two-dimensional signal g(,z, t) by time shifting the signals g,(t) for2 < i < n of said filtered 
ihift; 

>y performing a time shift on the signals g,(t) for 
re components of velocity Vtfc in the signals g,(t) for 1 < i < n ; 

b) filtering said shifted signal using the method according to any one of claims 1 to 6 ; 

c) producing a second two-dimensional signal g a (z, t) by time shifting the signals g,(t) for 2 < i < n of said 
filtered signal in order to cancel the preceding time shift ; and 

C) taking the difference between the signals g t (z, t) and g 2 (z, t) in order to produce a signal g, 2 (z, t) which 
is predominately constituted by waves having velocities lying between Vc, and Vc* 

8. Aseismic exploration installation for implementing the method according to claim 1 , the installation com- 
prising at least one sound wave detector together with means for displacing said detector along a borehole, 
means for recording a two-dimensional signal g(z, t) constituted by the signals gi(t) for 1 < i < n produced by 
the detector, said signals including upgoing wave components and downgoing wave components, and velocity 
filtering means for separating said components, the installation being characterized in that said velocity filtering 
means comprise : 

a memory (2 



6) for applying the Hilbert operator H t with respect to the variable t to the signal deli- 



calculating means (28) for determining the value of the coefficient B as a function of the signals delivered 
by the memory (22) and by the preceding calculating means (26) ; 

a multiplier (30) for multiplying together the coefficient s, the coefficient B, and the signal produced by the 



an adder (32) receiving the signal contained in the memory (22) and the signal produced by the multiplier 
(30); 

a divide-by-two circuit (34) receiving the signal produced by the adder (32) ; 

a multiplexer (38) having one input connected to the output from the divide-by-two circuit (34) and having 
two outputs, one of which is connected to a memory (40) and the other of which is connected to the cal- 
culating means (24) ; 

a counter (36) controlling the multiplexer, said counter controlling the number of times the operator A is 
applied to the signal g(z, t) ; 



g(z + 2m-l) - g(z-2m + l) 



le function and [ is an integer. 



1 o. An installation according to claim 9, characterized in that I is ec 
• ■ ad in that Msec 



12. An installation according to any of claims 8 to 1 1 , characterized in 
ments the operator H, defined by : 

HJg(t)] = IFF D.sgn(a.).FTfe(t)] 



IFTreL 

sgn(m) designates the sign of the variable a, which is 
| is the complex operator such that f = - 1 . 

13. An installation according to any one of claims 8 to 12, characterized in that the counter (36) cor 
the multiplexer (38) so that the operator A is applied only once to the signal g(z, t). 

14. An installation according to any one of claims 8 to 13, characterized in that it further includes ; m 
for time shifting the signals g,(t) for 2 < l< n constituting the signal g(z, t) stored in the memories (22, 
together with means for subtracting the signals g,(z. t) and g 2 (z. t) contained in the memory (40) from 



PatentansprQehe 

LEinseismisches 



ie es erzeugt wird durch mindestens einen seismischen Wel- 

z, z„ in einem Bohrioch plaziert ist, wobei jedes der 

" an und die Signale aufwarts 



A=1/2pd + e.B.H,.DJ 



in Wert abhfingt von dem Signal, auf das der Operator A angewandt 



ezuglich der Variable n t ist ; 



Dein Differenzierungsoperator bezuglich der Variablen z i=i , »,» 
wei^roperf^ 

Signal Oberwiegend gebildetwird durch aufwarts gehende Wellen, wenn s = + 1 ist, und durch abwartsgehende 
inzeichnet, daB der Operator D z definiert ist durch 



2. Ein Verfahren nach Anspruch 1 , d£ 



1 



• z 



n-1 

worin g_ eine Funktion ist und I eine ganze ZaW ist ' ' 
3 Ein Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB I gleich 1 ist. 

4. Ein Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB 1 gleich 2 list ' 

5. Ein Verfahren nach einem der Anspruche 1 bis 4. dadurch gekennzeichnet, daB der Operator H, definiert 
ist durch: 



HJg(t)] = I FT B-sgn((»).FT(g(r.)] 
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g eine Funktion ist ; 

FT die Fourier-Transformation reprasentiert ; 
IFT die inverse Fourier-Transfomiation re 
sgn(<g) das Vorzeichen der Variablen to bezeichnet, was das Duale von t ist ; und 
I der Komplexoperator ist derart, daR p = - 1 ist. 

mal auf das Signal g(z, t) angewandt wird. 

7. Ein Verfahren der Gesohwindigkeitsfilterung eines zweidimensionalen Signals g(z, t), aufgebaut aus ei- 
nem Satz von Signalen g,(t) fur 1 < i < n, welche Signaie erzeugt werden durch Schallwellendetektoren, die an 

unterschiedlichen Tiefen z, z z„ im Teufenbereich eines Bohriochs plaziert sind, wobei jedes Signal 

erzeugt wird im Ansprechen auf Sohallwellen, die von eineran der Erdoberflache befindlichen Quelle emittiert 
werden, wobei die Gesohwindigkeitsfilterung dazu dient, die Wellen in dem Signal g(z, t) zu verstarken mit 
Geschwindigkeiten, die zwischen vorbestimmten Gesohwindigkeiten Vc, und V02 liegen, welches Verfahren 
dadurch gekennzeichnet ist, daR es die folgenden Sohritte umfaRt : 

A) a) Erzeugen eines ersten verschobenen zweidimensionalen Signals durch Anwenden einer Zeitver- 

schiebung auf das Signal g,(t) fur 2 < i < n zwecks zeitlicher Ausfluchtung der Wellenkomponenten der 

Geschwindigkeit Vc, in den Signalen gfi) fur 1< i < n ; 

b) Filtem der verschobenen Signaie unter Verwendung des Verfahrens nach einem der Anspruche 1 bis 
6; 

c) Erzeugen eines ersten zweidimensionalen Signals g(,z, t) durch Zeitverschlebung der Signals g,(t) fur 



sohiebung auf die Signaie a(t) fur 2 < i < n zwi 
Geschwindigkeit VC2 in den Signalen g,(t) fOr 1 < I < n ; 

b) Filtern des verschobenen Signals unter Verwendung des Verfahrens nach einem der Anspruche 1 bis 

6; 

c) Erzeugen eines zweiten zweidimensionalen Signals g 2 (z, t) durch Zertverscliiebung der Signaie g,(t) fur 
2 < i < n des gefllterten Signals zwecks Ausioschung der vorhergehenden Zeitverschlebung, und 

C) Verwenden der Differenz zwischen den Signalen g,(z,t) und g 2 (z, t) zwecks Erzeugung eines Signals 
g 12 (z, t), das Oberwiegend von Wellen gebildet wird mit Geschwindigkeiten. die zwischen Vc, und Vcj lie- 
gen. 

8. Seismische Explorationsanlage fur die Durchfuhrung des Verfahrens nach Anspruch 1 , welche Anlage 
mindestens einen Schallwellendetektor mit Mitteln fur die Veriagerung des Detektors langs eines Bohriochs 
umfaRt, Mittel fur das Aufzeichnen eines zweidimensionalen Signals g(z, t) umfaRt, gebildet durch die Signaie 
g,(t) fur 1< i < n, erzeugt von dem Detektor, welche Signals aufwarts gehende Wellenkomponenten und 
abwarte gehende Wellenkomponenten umfassen, und Geschwindigkeitsfiltermittel fur das Trennen der 
genannten Komponenten, welche Anlage dadurch gekennzeichnet ist, daR die Geschwindigkeitsfiltermittel 
umfassen : 

einen Speicher (22) fur den Empfang des zweidimensionalen Signals g(z, t) ; 

Berechnungsmittel (24) fur die Anwendung des diskreten Differenzierungsoperators D z bezuglich der 
Variablen z in dem in dem Speicher (22) enthaltenen Signal ; 

Berechnungsmittel (26) fur das Anwenden des Hitbert-Operators H, bezuglich der Variablen t auf das 
Signal, geliefertvon den vorhergehenden Berechnungsmitteln (24) ; 

Berechnungsmittel (28) fur die Bestimmung des Wertes des Koeffizienten B als Funktion der Signaie, gelie- 
fertvon dem Speicher (22) und von dem vorhergehenden Berechnungsmittel (26) ; 
eine Multipliziereinrichtung (30) fur das Multiplizieren des Koeffizienten 5 des Koeffizienten B und des 
Signals, erzeugt durch die Berechnungsmittel (26) ; 

eine Addiereinrichtung (32), die das in dem Speicher (22) enthaltene Signal und das von der Mi 
einrichtung (3 




auf das Signal g(z, t) angewandt wird ; 1 



einen Speicher (40) ange- 
wie oft der Operator A 



D 2 [g<=o] - J2 



g(z + 2m-l) - sr(2-2m+l) 



worin g_ eine Funktion und I eine ganze Zahl ist 

10. Eine Anlage naoh Anspruch 9, dadurch gekennzeichnet, da& I gleich 1 ist. 

11. Eine Anlage naoh Anspruch 9. dadurch gekennzeichnet, daB \ gleich 2 ist. 

12. Eine Anlage nach einem der Anspruche 8 bis 11, dadurch gekennzeichnet, dad die Berechnungsmittel 
(26) den Operator H, implementieren, definiert durch : 

H 1 Ig(t)] = IFT[j.sgn( ( D).FT(g(t)] 

g eine Funktion ist; 

FT die Fourier-Transformation reprasentiert ; 

I FT die inverse Fourier-Transformation reprasentiert ; 

sgn(m) das Vorzeichen der Variablen <o bezeichnet, die das Dual von t ist ; und 

j der Komplexoperator ist derart, daB j 2 = - 1 . 

1 3. Eine Anlage nach einem der Anspruche 8 bis 12, dadurch gekennzeichnet, daB der ZShler (36) den 
Multiplexer (38) derart steuert, daB der Operator A nur einmal auf das Signal g(z, t) angewandt wlrd. 

14. Eine Anlage nach einem der Anspruche 8 bis 1 3, dadurch gekennzeichnet, daB sie ferner Mittel umfaBt 
fur die Zeitverschiebung der Signal g,(t) fur 2 < i < n, welche das Signal g(z, t), abgespeichert in den Speichem 
(22, 40), bildet, zusammen mil Mitteln fiir das Subtrahisren der Signale g,(z, t) und g 2 (z. t), enthalten im Spei- 
cher (40) vonelnander. 



1. Procede d'exploration sismique dans lequei 6n filtre un signal bidimensionnel g(z, t) compose d'un 
ensemble de signaux g,(t) oil 1< i< n, produits par au moins un detecteur d'ondes sismiques plac§ a diff§rentes 
profondeurs z,, ... z z„ dans un puits, chaque signal §tant produit en reponse a une emission d'ondes sis- 
miques, lesdits signaux comprenant des composantes d'onde remontante et des composantes d'onde descen- 
dante, ledit filtrage consistent a separer lesdites composantes et etant caracterise en ce qu'on applique au 
signal bidimensionnel g(z, t) I'operateur A defini par : 

A =1/2 0d + e.B.H,.DJ 

ou 

Id es t I'operateur identite, 

— B est un facteur de normalisation, dont la valeur depend du signal auquel I'operateur A est applique, 

— H, est I'operateur de Hilbert unidimensionnel relatif a la variable t, 

— D est un operateur de differentiation relatif a la variable, z, et 
— s est un entier egal a + 1 ou - 1, 

ledit operateur etant applique de maniere recursive audit signal bidimensionnel g(z, t), le signal resultant 
comprenantdefacon predorninahte les ondes remontantes lorsque e = + 1,et les ondes descendantes lorsque 
e = -1. 

2. Procede selon la revendication 1, caracterise en ce que I'operateur D z est defini par : 



EP 0290 336 B1 



(g Cz) ) =2, — 



HJg(t)] = IFTD-sgn(o>).FT(g(t)] 

oil g est une fonction queloonque 

FT represents la transformation de Fourier 
I FT represente la transformation de Fourier inverse 
sgn(m) note le signe de la variable o, duale de t 
j est le nombre complexe tel que j 2 = - 1 . 

6. Precede selon I'une queloonque des revendications 1 a 5, caracterise en ce que I'operateur A est appli- 
que une seule fois au signal g(z, t). 

7. Precede de filtrage de Vitesse d'un signal bidimensionnel g(z, t) compose d'un ensemble de signaux 

gi(t), oil 1 < 1< n produits par detecteurs d'ondes acoustiques places a differentes profondeurs z,, z , z n 

dans un puits, chaque signal etant produit en reponse a une emission d'ondes acoustiques a partir d'une source 
en surface, pour renforcer les ondes du signal g(z, t) ayant une Vitesse comprise entre des vitesses pr£dSter- 
minees Vc, et Vc* ledit precede etant caracterise en ce qu il comprend les operations suivantes : 

A)a) produire un premier signal bidimensionnel decale en effectuant un decalage temporal des signaux 
gi(t), 2 < i < n, pour aligner temporellement dans les signaux g,(t), 1< i < n, les composantes des ondes 
de Vitesse Vc„ 

is 1a 6, 



n, dudit signal filtre, pour annuler le decalage tempore! precedent, 



2<i<n, p. 
seVc 2 , 



b) filtrer ledit signal decale se 

c) produire un 
2 < i< n, dudi _ 

C) faire la difference des signaux g, (z, t) et g 2 (z, t) pour produire un signal g, 2 (z, t) comprenant de facon 
predominate les ondes ayant une vitesse comprise entre Vc, et Vcj. 

8. Installation d'exploration sismique pour la mise en oeuvre du precede selon la revendication 1 , compre- 
nant au moins un detecteur d'ondes acoustiques et des moyens pour deplacer ce detecteur dans un puits, des 
moyens pour enregistrer un signal bidimensionnel g(z, t) constitue par les signaux g,(t), 1 < i < n, produits par 
le detecteur, lesdits signaux comprenant des composantes d'onde descendante et des composantes d'onde 
remontante, et des moyens de filtrage en vitesse pour separer lesdites composantes, caracterise en ce que 
lesdits moyens de filtrage comprennent : 

— une memoire (22) pour recevoir le signal bidimensionnel g(z, t), 

— un moyen de calcul (24) pour appliquer I'operateur D z de differentiation discret relativement a la variable 
z au signal contenu dans la memoire (22), 

— un moyen de calcul (26) pour appliquer I'operateur H, de Hilbert relativement a la variable t au signal 
delivre par le moyen de calcul (24) precedent, 

— un moyen de. calcul (28) pour determiner la valeur du coefficient B en fonction des signaux delivres par 



— un multipiieur (30) pour multiplier entre eux le coefficient s, le coefficient B et le signal produit par le 
moyen de calcul (26), 

— un additionneur (32) recevant le signal contenu dans la memoire (22) et le signal produit par le multipiieur 



14 
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le nombre de fois oil I'operateur 
calcul (24) met en oeuvre I'ope- 



rateur D z defini par : 



— g<2+2m-1) - gCz-2m+1) 

•jiw) -z: — ^ 



j g est une fonction queloonque et I est un nombre entie 

10. Installation selon la revendication 9, caracterisee 

11. Installation selon la revendication 9, caracterisee en ce que I est egal a 2. 



dications 8 a 1 1 , caracterisee en ce que le moyen de calcul 



6) met en oeuvre i'operateur H, definl par : 

HIgWl = IFT B.sgn(o>).Fr(g(t)] 



, le signe de la variable oj, duale de t 

InstaMton s^n^e^elconque des revendications 8 a 12, caracterisee en ce que le compteur (36) 
multiplexer (38) pour que I'operateur A ne soit applique qu'une seule fois au signal g(z, t). 
.altan selon I'une queloonque des revendications 8 i 13 earache an ce qu'elle comprend en 
, constituantle signal g(z, t) memonse 
,t),g 2 (z,t)contenusdanB 
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